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Large-Eddy Simulations (LES) of two-phase turbulent flows exhibit 
quantitative differences in particle statistics if compared to Direct Nu- 
merical Simulations (DNS) which, in the context of the present study, 
is considered the exact reference case. Differences are primarily due 
to filtering, a fundamental intrinsic feature of LES. Filtering the fluid 
velocity field yields approximate computation of the forces acting on 
particles and, in turn, trajectories that are inaccurate when compared 
to those of DNS. In this paper, we focus precisely on the filtering er- 
ror for which we quantify a lower bound. To this aim, we use a DNS 
database of inertial particle dispersion in turbulent channel flow and we 
perform a-priori tests in which the error purely due to filtering is singled 
out removing error accumulation effects, which would otherwise lead to 
progressive divergence between DNS and LES particle trajectories. By 
applying filters of different type and width at varying particle inertia, 
we characterize the statistical properties of the filtering error as a func- 
tion of the wall distance. Results show that filtering error is stochastic 
and has a non-Gaussian distribution. In addition, the distribution of 
the filtering error depends strongly on the wall-normal coordinate being 
maximum in the buffer region. Our findings provide insight on the ef- 
fect of subgrid-scale velocity field on the force driving the particles, and 
establish the requirements which a LES model must satisfy to predict 
correctly the velocity and the trajectory of inertial particles. 
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I. INTRODUCTION 



The dispersion of small inertial particles in inhomogeneous turbulent flow 
has been long recognized as crucial in a number of industrial applications and 
environmental phenomena. Examples include mixing, combustion, depulver- 
ation, spray dynamics, pollutant dispersion or cloud dynamics. In all these 
problems accurate predictions are important, yet not trivial to obtain because 
of the complex phenomenology that controls turbulent particle dynamics. 

Direct Numerical Simulations (DNS) of turbulence coupled with Lagrangian 
Particle Tracking (LPT)^"- have demonstrated their capability to capture the 
physics of particle dynamics in relation with turbulence dynamics and have 
highlighted the key role played by inertial clustering and preferential concen- 
tration in determining the rates of particle interaction, settling, deposition 
and entrainment. Due to the computational requirements of DNS, however, 
analysis of applied problems characterized by complex geometries and high 
Reynolds numbers demands alternative approaches, among which Large-Eddy 
Simulation (LES) is gaining in popularity, especially for cases where the large 
flow scales control particle motion (e.g.-*^). 

LES is based on a filtering approach in which the unresolved Sub-Grid 
Scales (SGS) of turbulence are modeled. A major issue associated with LES 
of particle-laden flows is modeling the effect of these scales on particle dynam- 
ics. Flow anisotropy adds a further complicacy. In particular, recent studies 
on particle dispersion in wall-bounded turbulenceii"— have demonstrated that 
use of the filtered fluid velocity in the equation of particle motion with no 
model for the effect of the SGS fluid velocity fluctuations leads to significant 
underestimation of preferential concentration and, consequently, to weaker de- 
position fluxes and lower near-wall accumulation. Even though several studies 
have demonstrated that predictions may be improved using techniques like 
filter inversion or approximate deconvolutioniiii^ii^"— , the amount of velocity 
fluctuations that can be reintroduced in the particle equations does not ensure 
a quantitative replica of DNS results. A further effort is required to model 
the scales smaller than the filter width but extremely significant for particle 
dynamics, which are inevitably lost in the filtering procedure^. Attempts to 
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model SGS effects on particles have been made using fractal interpolationr^ 
or kinematic simulations^^ in channel flow. Also these approaches, however, 
do not ensure complete removal of inaccuracies in the quantitative prediction 
of local particle segregation and accumulation (particularly in the near-wall 
region). Stochastic models have also been proposed, yet we are only aware of 
validations in homogeneous isotropic turbulence^>^>^: these models have still 
to be assessed in non-homogeneous anisotropic flows. 

To devise accurate and reliable SGS models for LPT, precise quantification 
of filtering effects on particle dynamics is clearly a necessary step. Aim of 
this paper is to propose a simple procedure to characterize the subgrid error— 
purely due to filtering of the fluid velocity field. This error, referred to as pure 
filtering error hereinafter, can not be avoided but potentially can be corrected. 
We thus focus on an ideal situation in which LES provides the exact dynamics 
of the resolved velocity field {ideal LES). In ideal LES, modeling and numer- 
ical errors (incurred because a real LES provides only an approximation of 
the filtered velocity) and interpolation errors (on coarse-grained domains, for 
instance )ii are assumed negligible, and a lower bound for the pure filtering 
error can be identified. Peculiar feature of the procedure proposed here is that 
it removes time accumulation of pure filtering errors on particle trajectories. 
Error accumulation originates from inaccurate estimation of the forces acting 
on particles obtained when the filtered fluid velocity is supplied to the equation 
of particle motion. As a result, particle trajectories in LES fields progressively 
diverge from particle trajectories in DNS fields, considered as the exact refer- 
ence for the present study: the flow fields seen by the particles become less and 
less correlated and the forces acting on particles are evaluated at increasingly 
different locations. We remark that the effects due to these two errors can not 
be singled out easily when comparing statistics of particle velocity and concen- 
tration obtained from LES to the reference DNS data. To remove this effect, 
we perform here a-priori tests in which: (i) particle trajectories are computed 
from the DNS fields, (ii) the DNS fields are coarse-grained through filtering 
and, (in) particles are forced to evolve in the filtered DNS fields along the 
DNS trajectories. A similar a-priori analysis has been performed recently to 
evaluate pure filtering error and time accumulation error for the case of tracer 
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particle dispersion in homogeneous isotropic turbulence^. In this paper, we 
focus on inertial particles in turbulent channel flow. 

A-priori evaluation of the statistical properties of the filtering error provides 
useful information about the key features that should be incorporated in SGS 
models for LPT to compensate for such error. Statistics are computed ap- 
plying filters of different type and different widths, corresponding roughly to 
varying amounts of resolved flow energy, and considering point particles with 
different inertia that obey one-way coupled Stokes dynamics. One point at 
issue is indeed whether SGS models for LPT should take into account inertial 
effects explicitly. Another point to be explored is whether the characteristics 
of SGS models for LPT should adapt "dynamically" to respond to the local 
inhomogeneity and anisotropy of the flow. Our study aims at addressing these 
issues through statistical characterization of the effects of the unresolved veloc- 
ity field on the force driving the particles. Possible ways to improve prediction 
of LES applied to turbulent dispersed flow should be benchmarked against this 
type of statistical information at the sub-grid level. 

The paper is organized as follows. The physical problem and the numerical 
methodology are described in Sec. [Tll The statistical moments and the proba- 
bility distribution functions of the filtering error are presented in Sec. IIIII and 
in Sec. llVt respectively. Finally, concluding remarks are given in Sec. |Vl 

II. PHYSICAL PROBLEM AND NUMERICAL 
METHODOLOGY 

A. Particle-laden turbulent channel flow 

The reference flow configuration consists of two infinite fiat parallel walls: 
the origin of the coordinate system is located at the center of the channel 
and the x—, y— and z— axes point in the streamwise, spanwise and wall- 
normal directions respectively. Periodic boundary conditions are imposed on 
the fluid velocity field in x and y, no-slip boundary conditions are imposed at 
the walls. In the present study, we consider air (assumed to be incompressible 
and Newtonian) with density p = 1.3 kg and kinematic viscosity u = 
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15.7x10"^ s^^. The flow is driven by a mean pressure gradient, and the 
shear Reynolds number is Rcr = Urh/u = 150, based on the shear (or friction) 
velocity, u^-, and on the half channel height, h. The shear velocity is defined 
as Ur = {tw/pY^"^, where is the mean shear stress at the wall. 

Particles are modeled as pointwise, non-rotating rigid spheres with density 
pp = 10^ kg m~^, and are injected into the flow at concentration low enough 
to consider dilute system conditions (no inter-particle collisions) and one-way 
coupling between the two phases (no turbulence modulation by particles). 
The motion of particles is described by a set of ordinary differential equations 
for particle velocity and position. For particles much heavier than the fluid 
{Pp/p ^ 1) Elghobashi and Truesdell^ have shown that the most significant 
forces are Stokes drag and buoyancy and that Basset force can be neglected 
being an order of magnitude smaller. In the present simulations, the aim is to 
minimize the number of degrees of freedom by keeping the simulation setting 
as simplified as possible; thus the effect of gravity has also been neglected. 
With these assumptions, a simplified version of the Basset-Bousinnesq-Oseen 
equation^^ is obtained. In vector form: 

dvp _ u^j22ii(i + 0.15iJe»^«»') . (2) 

CIl 

In Eqns. ([I]) and ([2]), Xp is particle position, Vp is particle velocity, = 
u(xp(t),t) is the fluid velocity at the particle position. Tp = Ppd^/lSp is the 
particle relaxation time, and Rcp = dp\\p — Us\p/p is the particle Reynolds 
number, dp and p being the particle diameter and the fluid dynamic viscosity 
respectively. 



B. DNS/LPT methodology 

The Eulerian flow field is obtained using a pseudo-spectral DNS flow solver 
which discretizes the governing equations (Continuity and Navier-Stokes for in- 
compressible flow) by transforming the field variables into wavenumber space, 
using Fourier representations for the periodic streamwise and spanwise direc- 
tions and a Chebyshev representation for the wall-normal (non-homogeneous) 
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direction. A two level, explicit Adams-Basliforth scheme for the non-linear 
terms, and an implicit Crank-Nicolson method for the viscous terms are em- 
ployed for time advancement. Further details of the method can be found in 
previous articles (e.g. Pan and Banerjee^^). Calculations were performed on 
a computational domain of size x Ly x Lz = Anh x 2'Kh x 2h in x, y and 
z respectively, corresponding to 1885 x 942 x 300 in wall units. Wall units 
are obtained combining Ur, v and p. The computational domain is discretized 
in physical space with A^^. x Ny x = 128 x 128 x 129 grid points (corre- 
sponding to 128 X 128 Fourier modes and to 129 Chebyshev coefficients in the 
wavenumber space). This is the minimum number of grid points required in 
each direction to ensure that the grid spacing is always smaller than the small- 
est flow scale and that the limitations imposed by the point-particle approach 
are satisfied.— 

To calculate particle trajectories a Lagrangian tracking routine is coupled 
to the fiow solver. The routine solves for Eqns. ([T]) and ([2]) using a 4*^-order 
Runge-Kutta scheme for time advancement and 6*^-order Lagrangian polyno- 
mials for fiuid velocity interpolation at the particle location. At the beginning 
of the simulation, particles are distributed randomly within the computational 
domain and their initial velocity is set equal to that of the fiuid at the particle 
initial position. Periodic boundary conditions are imposed on particles moving 
outside the computational domain in the homogeneous directions, perfectly- 
elastic collisions at the smooth walls are assumed when the particle center is 
at a distance lower than one particle radius from the wall. For the simula- 
tions presented here, large samples of 10^ particles, characterized by different 
response times, are considered. When the particle response time is made di- 
mensionless using wall variables, the Stokes number for each particle set is 
obtained as St = Tp/rj where Tf = v/u^ is the viscous timescale of the fiow. 
Three different sets of particles corresponding to values of the Stokes number 
St = 1, 5 and 25 have been considered in this study: Table [T] shows the relevant 
physical parameters of each set. The total tracking time in dimensionless wall 
units (identified hereinafter by superscript +) is AT'^ = 21150 for all particle 
sets^, the timestep size for particle tracking being equal to the timestep size 
used for the fiuid: 5t+ = 0.045. 
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C. A-priori LES methodology and computation of the filtering 
error 



In the a-priori tests, LPT is carried out replacing in Eq. ([2]) with the 
filtered fluid velocity field, u(xp,t). This field is obtained through explicit 
filtering of the DNS velocity by either a cut-off filter or a top-hat filter. Both 
filters are applied in the wave number space to velocity components in the 
homogeneous streamwise and spanwise directions: 

\ G{ki) ■ G{k2) ■ Ui{Ki,K2,z,t) if < |kc| with j = 1,2 
Ui{x,t) = FT < (3) 

I otherwise 

where FT is the 2D Fourier Transform, = tt/A is the cutoff wave number 
(A being the filter width in the physical space), Ui^ni, k,2i z,t) is the Fourier 
transform of the DNS fluid velocity field, namely Ui{K,i, k,2, z, t) = -FT[Mj(x, t)] 
and G{K,j) is the filter transfer function: 

{1 for the cut-off filter 

sin(KA/2) 
top-hat filter 

Three different filter widths are considered, which provide a grid Coarsening 
Factor (CF) in each homogeneous direction of 2, 4 and 8 with respect to DNS, 
corresponding to 64 x 64, 32 x 32 and 16 x 16 Fourier modes respectively. 
Note that CF=2 and CF=4 yield grid resolutions that are commonly used in 
LES, whereas CF=8 corresponds to a very coarse resolution characteristic of 
under- resolved LES. Data are not filtered in the wall-normal direction, since 
the wall-normal resolution in LES is often DNS-like.— 

The pure filtering error is computed under the ideal assumption that all 
further sources of error affecting particle tracking in LES can be disregarded 
and that the exact dynamics of the resolved velocity field are available. Errors 
due to the SGS model for the fluid and to numerics are thus neglected and a 
lower bound for the filtering error is identified by removing error accumulation 
effects due to progressive divergence between DNS and LES particle trajecto- 
ries. Time accumulation of the pure filtering error is prevented by computing 
particle trajectories in DNS fields and then forcing particles to evolve in fil- 
tered DNS fields along the DNS trajectories. At each time step n and for each 
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particle k we thus impose: 

where superscripts LES and DNS identify particle position and velocity in 
a-priori LES and in DNS, respectively. The pure filtering error made on the 
k-th particle at time is computed as: 

6u = 6u{^p,k{n,n = u(xj,,fc(r),r) - u(xp,fc(r),r) = u, - a, . (6) 

Due to the wall-normal dependence of the fluid velocity (and of in partic- 
ular), it is not straightforward to select a characteristic velocity to normalize 
6u. We thus decided to provide measures of the absolute filtering error rather 
than the percent filtering error, using Eq. (I6l) to quantify exactly the filtering 
error affecting the fluid velocity seen by the particles. This quantification relies 
on an idealized situation which serves the purpose of isolating the contribution 
due to filtering of the velocity field to the total error associated with LPT in 
LES. Note also that 6u is computed along the "exact" trajectory of each parti- 
cle, that is at x^^'^(t"'). As shown in the following, the corresponding behavior 
observed for 6u may be significantly different from the Eulerian measure of the 
filtering error usually found in the literature, which is obtained as difference 
between DNS and LES velocities at fixed points in space. 

Computation and statistical characterization of 6u were carried out for all 
particle sets reported in Table [land applying filters of different type and width. 
The averaging time window required for converged statistics was At~^ ~ 8285. 



III. STATISTICAL MOMENTS OF THE PURE FILTERING 
ERROR 

The Eulerian statistical moments of the filtering error, quantified here by 
its lower bound 6u (Eq. ([6])), are presented and analyzed in this section. 
Statistics were computed by dividing the computational domain in wall-parallel 
slabs having dimensions L^, Ly and Az{i) with i = l,Nz, where Az{i) is 
the difference between two adjacent Chebyshev collocation points. Ensemble 
averaging is carried out over particles that, at a given time instant, are located 
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inside the same slab; time averaging is performed over the time interval At^ . 
All quantities shown in this section are non dimensional and expressed in wall 
units. 

A. Influence of particle inertia on the filtering error 

To analyze the effect of particle inertia on the behavior of 6u, we compare 
statistics relative to all particle sets but we focus one filter type (cut-off) and 
one filter width (CF=4). This filter width is chosen as it removes a significant 
amount of spatial information on the turbulent structures which interact with 
particles without generating overly under-resolved fields. The effects due to 
different filter types and widths will be investigated in Sec. IIIIBI 

Figure [T] shows the mean absolute value {6u) of the filtering error along 
the streamwise, spanwise and normal directions as a function of the distance 
from the wall. Hereinafter (■) identifies time- and ensemble-averaged quanti- 
ties. The qualitative behavior of the streamwise component (Sux) (Fig. [1^) is 
similar for all particle sets: profiles appear rather flat in the central region of 
the channel (60 < < 150) where {Su^) is nearly zero and develop a peak 
in the near-wall region where {6ux) becomes more and more negative as the 
wall is approached. Negative {6ux) means over-estimation of the streamwise 
drag acting on the particles due to filtering. Maximum overestimation occurs 
at ^ 15 — 20: the location of this maximum does not change for different 
inertia. Inertia has a strong effect on the magnitude of the peak, which is max- 
imum for St = 5 and falls off on either side of this value. This non-monotonic 
dependence can be explained considering that inertial particles are low-pass 
filters that respond selectively to removal of sub-grid flow scales according to 
a characteristic frequency proportional to l/r^r^. When the frequency of the 
removed scales is equal to the characteristic frequency of the particles, particle 
response is maximized. This is precisely the case of Fig. [1^), in which filtering 
with CF=4 are observed to produce larger errors for the St = 5 particles com- 
pared to the other particle sets. We remark that inertia also plays a role in 
transferring filtering errors on fluid velocity into filtering errors on drag, which 
are proportional to (5u) ■ St. 
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Filtering has no significant effect in mean on the spanwise component {5uy) 
(Fig. [T|d) which is characterized by very small values — (9(10~^) —oscillating 
around zero along the entire channel height. In the spanwise direction, how- 
ever, both (us) and {us) attain negligible values. More interesting observations 
can be drawn examining the wall-normal component, {5uz) (Fig. [ifc). First, 
we notice that {5uz) profiles develop a positive peak in the near-wall region 
which increases in magnitude with the Stokes number. Moving away from this 
region, values become negative and then relax to zero toward the center of the 
channel. This means that drag may be either underestimated or overestimated 
depending on the particle-to-wall distance, a feature not trivial to be repro- 
duced in a model. The observed behavior of {5uz) also shows that filtering in 
the homogeneous directions has consequences on the unfiltered inhomogeneous 
direction. The filtering error propagates in the wall-normal direction according 
to the relation ^^(x, t) = FT"-*^ G{k,i) ■ G{k,2) ■ Uz^t^i, ti2^ z,t) , see Eq. ([3]). 

We remark that all components of 5\i would be rigorously equal to zero if 
computed at fixed grid points rather than along particle trajectories. This 
because filtering does not affect the mean value of the velocity. In a-priori and 
a-posteriori tests, it is customary to analyze the effects of filtering on the LES 
fluid velocity fields at fixed points (see e.g.— When dealing with LPT, 
however, a more natural approach is to evaluate filtering effects on the fluid 
velocity seen by the particles: in this way, both filtering effects on turbulent 
structures and preferential concentration effects are reproduced in the behav- 
ior of 5\i. From a physical viewpoint, this establishes a direct link between 
5\i and the phenomenology of particle dispersion in turbulent boundary layer. 
To elaborate, it is well known that near-wall turbulence is characterized by 
the presence of low-speed and high-speed streaks and that inertial particles 
preferentially sample low-speed streaks^'^'^. Visual evidence of this tendency 
is provided in Fig. |2^, in which the instantaneous position of the St = 25 
particles comprised within a near-wall region 10 wall units thick is superposed 
to fluid streaks on a (x; t/)-plane close to the wall {z~^ = 10). Fluid streaks are 
rendered using colored contours of the streamwise fluid velocity fluctuation u'^ 
in a DNS flow field. The effect of filtering becomes clear upon comparison of 
Fig. [2^ with Figs. [2]D-[2t, which demonstrate how rendering of the streaky 
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structures deteriorates when using a coarse-grained flow field (with coarsening 
factors 4 and 8, in this specific case). Filtering attenuates the velocity fluctu- 
ations and smooths the streaks thus increasing the fluid velocity seen by the 
particles: this is precisely what causes the negative peaks of {5ux) in Fig. [1^. 
The behavior of {5uz) can also be explained in terms of filtering effects on the 
wall-normal fluid velocity seen by the particles, shown in Fig. |3] for all parti- 
cle sets. Lines refer to DNS profiles {{ug^z)), symbols to a-priori LES profiles 
{{us,z))- Note that (Su^) = {us^z — Us,z) = (^s,z) ^ {us,z)- Even though parti- 
cles exhibit a mean wallward drift, both (us^z) and (us^z) are positive (namely 
directed away from the wall) in the near-wall region due to particle preferen- 
tial concentration in ejection-like environments, where low-momentum fluid is 
entrained off the wall^. Filtering reduces the wall-normal fluid velocity seen 
by the particles very close to the wall while increasing it towards the center 
of the channel. This is due to smoothing of the near-wall structures and at- 
tenuation of the sweep events, which are characterized by coherent motions of 
high-momentum fluid to the wall (see e.g. Fig. 6 in Marchioli et al.— ). 

Examination of {6u) confirms that LES is inaccurate in capturing the near- 
wall turbulent structures and in reproducing quantitatively their interaction 
with particles. Errors introduced by filtering lead to underestimation of pref- 
erential concentration and wall accumulation, as demonstrated by previous 
studiesii"— . To provide a thorough characterization of these errors, in Figs, 
m and [5] we examine the higher order moments. Analysis of the root mean 
square (rms) of the filtering error, shown in Fig. HI highlights the reduction of 
the fluid velocity fluctuations in LES already found in error analyses carried 
out at fixed (Eulerian) points^. SGS models previously used for LPT in LES 
fields, like approximate deconvolution or filtering inversion, are mainly aimed 
at counteracting this effect by recovering the correct amount of fluid and par- 
ticle velocity fluctuations in the particle equationsi^ii^. Note that, except in 
the very-near wall region, rms values of the error are in general much higher 
than the mean values along the same flow direction: only in the streamwise 
direction and within the buffer layer, rms{ux) and {5ux) have the same order 
of magnitude. No significant effects due to particle inertia are observed. 

The skewness and flatness factors of the filtering error in each flow direction 
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are shown in Fig. [5l indicated as ^(^u) and F{6u.) respectively. As for the rms, 
the most interesting results are obtained for the skewness in the streamwise 
and wall-normal directions {S{6ux) and S{6uz) in Figs. |5^ and |5^, respec- 
tively): values of the skewness in the spanwise direction {S{6uy) in Figs. [5t 
remain very small along the entire channel height, roughly oscillating around 
zero. Qualitatively, the skewness profiles in the streamwise and wall-normal 
directions are similar: moving from the wall toward the center of the channel, 
both S{6ux) and S{6uz) exhibit high positive values in the viscous sublayer 
followed by low negative values in the buffer region which then relax toward 
zero at the channel centerline. This behavior indicates that significant asym- 
metric deviations from Gaussianity are expected in the probability distribution 
functions of the streamwise and wall-normal components of the filtering error 
due to significant contributions given by strong isolated events. This finding 
is consistent with previous analyses^iii^ which showed that filtering impacts 
particle dynamics mainly through attenuation of low-speed streaks (negative 
velocity fiuctuations) and ejection events. Regarding inertial effects, changes 
in the Stokes number modify the S{6ux) profiles particularly in the central 
portion of the channel, while larger deviations among the S{6uz) profiles are 
observed in the near-wall region. Inertia has almost no effect on the flatness 
factors, shown in Figs. [5]d, |5]i), and[5]F: all proflles are characterized by high 
values at the wall which gradually decrease to centerline values of about 5 for 
F{6ux) and F{6uz) and of about 4 for F{6uy). Note that these values are 
higher than 3, which would be obtained if the error distribution was Gaussian. 

B. Influence of filter type and width on the filtering error 

To investigate the effect of fllter type and width on the statistical moments 
of 6u, we focus on the St = 5 particles. Filter effects on the mean value 
of the flltering error {6u) are shown in Fig. [61 The streamwise component 
{{6ux), Fig. [6^), exhibits a near- wall negative peak whose absolute value 
increases with the fllter width, namely with CF. This behavior is common 
to all particle sets investigated in this study. Regarding the fllter type, it 
is observed that the top-hat fllter introduces a slightly stronger smoothing 
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of the fluid velocity fluctuations compared to the cut-off filter, especially for 
CF=2 and CF=4, The mean spanwise component of {6u) is not shown since 
it is characterized by values that are nearly zero throughout the channel. The 
wall-normal component (Suz) is shown in Figure [6]3. It is evident that the 
behavior of (Suz) is influenced significantly by the filter width rather than by 
the filter type. As CF increases, the profiles develop positive and negative 
peaks of different intensity; also the distance from the wall at which these 
peaks occur changes considerably: all peak locations shift toward the center of 
the channel, indicating that filtering errors propagate their effect over a larger 
portion of the computational domain. 

Analysis of the higher-order moments also shows that effects due to the filter 
width are much stronger than those associated with the filter type. No effect 
of the filter type is observed on the rms of the filtering error, reported in Fig. 
[71 Profiles change only with the filter width, reaching higher values as CF 
increases because of the larger amount of fluid velocity fluctuations damped in 
more coarse-grained fields. Filter type effects are minor also on the skewness 
and flatness factors, shown in Fig. [HI When the filter width is increased, the 
skewness factor in the streamwise direction (Fig. [HH) increases to high positive 
values very close to the wall and decreases to low negative values outside the 
buffer layer {z'^ > 30). In this outer region of the flow, the skewness in the 
wall-normal direction (Fig. [8^) increases significantly, switching from negative 
to positive values as CF increases. High positive value of the skewness mean 
that 6u more frequently attains large positive rather than negative values, and 
of course the reverse for negative skewness. The largest filter width effects 
on the flatness factor are observed in the streamwise direction (Fig. [8b) and 
produce a decrease of F{6ux) throughout the entire channel: for small and 
intermediate widths {CF = 2 and 4) the flatness factor is everywhere higher 
than 3, for large widths (e.g. CF = 8) values fall below 3 in the region 
10 < 2;"'" < 30. We remind that high values of the flatness factor, like those 
characterizing the spanwise direction (Fig. [8)1), and the wall-normal direction 
(Fig. [Hf), indicate that error fluctuations are often larger than the variance of 
the distribution and that error fluctuations have an intermittent character. 
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IV. PROBABILITY DISTRIBUTION OF THE FILTERING 
ERROR 

In the previous section, we have characterized statistically a lower bound for 
the error produced by the filtering procedure inherent to the LES approach. 
This error is intrinsic of LPT in LES flow fields and its mean and higher order 
moments depend both on particle inertia and on filter width. In this section, 
we focus on the characterization of the one-point Eulerian probability density 
function (PDF) of the filtering error. This function can be derived from the 
more general Lagrangian PDF and is the most practical for direct evaluation 
of the errors made in the calculation of the statistical moments^^i^E. In the 
frame of the present study, the issues of interest about the PDF of the filtering 
error can be summarized by the following questions: 

1. Is the filtering error deterministic or stochastic? 

2. Does it have a simple distribution, for instance Gaussian stochastic? 

3. Is it universal or flow-dependent ? 

To answer these questions, we consider the PDFs of the streamwise and wall- 
normal component of the lower-bound filtering error 6u. These quantities, 
indicated hereinafter as p{6ux) and p{6uz) respectively, have been computed 
averaging over the same fluid slabs defined in Sec. IIIII and over the same 
time interval given in Sec. Ill CI In addition, PDFs have been normalized 
and shifted to be centered around zero: without shifting the profiles, both 
p{Sux) and p{Suz) would be centered around different values of 6ux and 6uz for 
the different particle sets, as a result of the dependence of the filtering error 
on particle inertia (see Sec. IIII Ap . Probability distributions in the spanwise 
direction are not shown as they are very similar to those obtained for p{6uz) 
and add very little to the discussion. 

We examine first the PDF behavior at varying filter widths, shown in Figs. 
[9] and [lOl To this aim, we focus on results obtained for the St = 5 particles 
when the cut-off filter is used (as for the statistical moments of 6u, effects on 
the PDF behavior due to the filter type appear negligible and hence will not 
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be discussed). In these figures profiles relative to four different distances from 
the wall are considered to emphasize flow anisotropy effects: clearly, choosing 
a non-homogeneous anisotropic flow as reference configuration is of particu- 
lar importance to devise general conclusions. For comparison purposes, the 
normalized Gaussian distributions (lines) with variance equal to that of the 
computed PDFs (symbols) at the same wall distance are also shown. A first 
(expected) result is that the PDF variance increases with the filter width in all 
flow directions and at all wall-normal locations examined. Another important 
finding is that, in all cases examined, the error is not deterministic, displays 
very complex features and covers a wide range of values (typical of probabilis- 
tic process) in all cases. Specifically, the PDFs are not Gaussian and their 
behavior varies strongly with the distance from the wall. This observation is 
not obvious since the error behavior might be quite different and even perfectly 
Gaussian in spite of the non-homogeneous nature of the flow. Figs. [9] and [TO] 
show that both p(5u^) and p{Suz) are characterized by larger flatness than 
the corresponding Gaussian distributions for all filter widths and at almost all 
wall-normal locations. Furthermore, both PDFs are skewed near the wall —at 
~ 4, Figs. [9^ and [TOk —and outside the buffer region —at ^ 65, Figs. 
[9ti and [TOl d. This is again related to poor description of coherent structures, 
e.g. low-speed streaks and turbulent vortices. At intermediate wall-normal 
locations, PDFs exhibit a nearly Gaussian behavior only at 2;"'' ^ 17 for all 
filter widths but the largest —Figs. [9b and [TOb : at z'^ ^ 37, where Reynolds 
stresses reach a maximum, the non-Gaussian behavior of the PDF profiles in- 
creases again —Figs. [Hb and[TOb. It can also be observed that (i) p{Suz) profiles 
become flatter at increasing filter widths (Fig. [TOl) . and (ii) filter width effects 
lead to higher skewness and mildly decreasing flatness of the PDF, consistently 
with results of Figs. [8^ and[8)D. For the largest filter width (CF=8), there is a 
change from positive skewness in the very near wall region {z~^ ^ 4) to negative 
skewness away from the wall. 

Additional information on p{6ux) and p{Suz) can be obtained examining 
how these quantities are modified by particle inertia. To this aim, we focus 
on results obtained using a cut-off filter with CF=4, shown in Figs. [TT] and 
[T2I A first observation is that the effect of inertia on the shape of the PDFs is 
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almost negligible, as we have tried to emphasize by centering all profiles around 
5ui = 0, where i = x,z. A general feature of the PDFs reported in Figs. [TTland 
[T2] is that they exhibit higher tails compared to the corresponding Gaussian 
distributions, even near the channel center. This is in line with the analysis 
of the fiatness profiles carried out in Sec. IIII A[ The fiatness is larger for 
the wall-normal probability distributions and varies mildly with the distance 
from the wall, the only exception being the viscous sublayer (see profiles at 
~ 4). In terms of relative importance of the PDF tails, p{Su^) profiles 
exhibit asymmetries towards negative values of 6ux which become particularly 
evident outside the viscous sublayer {z~^ ^ 17, 37 and 65); whereas p{Suz) 
profiles are all symmetric with just a slight asymmetry observed at ~ 17. 
These findings are in agreement with the skewness profiles shown in Figure El 



Some general conclusions can be drawn from collective analysis of Figs. [9] 
to [121 First, intrinsic filtering errors of LES in LPT can be described as a 
stochastic process almost everywhere non-Gaussian with important and long 
tails (namely large fiatness) and appreciable skewness. This process is signif- 
icantly infiuenced by fiow anisotropy as demonstrated by the dependency of 
skewness of p{Sux) on the distance from the wall. This result is important 
from the modeling viewpoint, as it implies that neither deterministic models 
nor stochastic homogeneous models have the capability to correct fully in- 
trinsic error made in the LES approach due to filtering. Second, the model 
should adapt dynamically over space to account for local inhomogeneity and 
anisotropy of the fiow and to reproduce the correct crossing trajectory ef- 
fect: this latter effect, in particular, is due to the existence of a mean relative 
particle-to-fiuid velocity (rather than an instantaneous one) associated with 
the lack of correlation between the fluid velocity seen by the inertial parti- 
cles and the velocity of fluid particles^S. Finally, results demonstrate that a 
mandatory requirement to scrutinize model performance is the capability to 
reproduce correctly all spatial turbulent structures (or at least their statistical 
effect). 
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V. CONCLUDING REMARKS 



In this work, we proposed a simple procedure to characterize the filtering 
error intrinsic of Lagrangian particle tracking in LES flow fields. Specifically, 
we focused on the error purely due to filtering of the fluid velocity fields seen 
by the particles, for which we quantified a lower bound. This "minimum" error 
can not be avoided but potentially can be corrected. To this aim, we consid- 
ered an ideal situation in which the exact dynamics of the resolved velocity 
field is available and time accumulation of the pure filtering error on particle 
trajectories is eliminated a-priori. Evaluation of the statistical properties of 
the filtering error has been performed applying filters of different type and 
different widths, corresponding roughly to varying amounts of resolved flow 
energy, and considering heavy particles with different inertia. In particular, 
mean and higher-order moments as well as one-point probability distribution 
functions have been analyzed to provide information about the key features 
required in SGS models for Lagrangian particle tracking to compensate for 
such error. 

The first interesting result is that the mean value of the filtering error com- 
ponents along particle trajectories is not zero, as it would be obtained per- 
forming an Eulerian analysis of filtering effects at fixed points. Non-zero mean 
values result from filtering effects on the coherent turbulent structures which 
characterize the near-wall region and are preferentially sampled by the inertial 
particles. Smoothing of the low-speed streaks due to filtering is associated 
with a near-wall negative peak of the mean filtering error in the streamwise 
direction. Filtering effects on the near-wall vortical structures, and in turn 
on sweep and ejection events, lead to even more complicated mean profiles in 
the wall-normal direction: in this case, a near-wall positive peak followed by a 
negative peak farther from the wall is observed. In both directions, particle in- 
ertia and filter width change significantly the qualitative behavior of the error. 
Mean value analysis confirms that the error introduced by filtering on inertial 
particle motion is not just yield by the reduction of fluid velocity fluctuations, 
a well-known effect pointed out in Eulerian error analyses in LES. This sup- 
ports the conclusion, drawn in our previous works^*^, that SGS model for the 
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particle equations can not ensure accurate prediction of particle preferential 
concentration and near-wall accumulation by only reintroducing the correct 
amount of fluid velocity fluctuations. Improved predictions of LES applied 
to turbulent dispersed flow necessitate suitable benchmarks against statistical 
information at the sub-grid level, and the present study establishes precisely 
the minimum benchmarking requirements. 

Joint analysis of the the higher-order statistical moments and of the prob- 
ability distribution functions shows that the filtering error distribution differs 
significantly from a normalized Gaussian distribution with the same variance. 
All error components are characterized by a flatness generally larger than the 
Gaussian value of 3. The streamwise and wall-normal components are also 
characterized by skewed distributions, not observed in the spanwise direction: 
both positive and negative skewness values can be attained depending on the 
distance from the wall. This indicates that error behavior is intermittent and 
strongly affected by flow anisotropy. Changes in the probability distributions 
are also found at varying filter width. Besides the expected increase of the vari- 
ance of the PDF when the filter become coarser, noticeable effects are observed 
on the skewness and flatness of the streamwise PDF and on the skewness of the 
wall-normal PDF. No significant effects are observed when either filter type or 
particle inertia are changed. 

Present results emphasize the stochastic nature of the pure filtering error. 
In our opinion, a possible way to correct this error could be represented by a 
generalized stochastic Lagrangian model^i^ in which the statistical signature 
of coherent turbulent structures (e.g. low-speed streaks, sweeps and ejections) 
on particle dynamics is accounted for explicitly. In this perspective, the fil- 
tering error characterization and quantification provided by the present study 
represents a valuable tool to understand how these structures are affected by 
filtering and how sub-grid scale fluid velocities impact the force driving the 
particles. Models of this type have already been developed in the context 
of RANS methods22^ followed by extensions to heuristic particle deposition 
modeling^ and by first attempts at explicitly modeling the statistical effects 
of structures on particle statistics in near-wall regions^. Much work remains 
to be done. The first necessary step is of course to extend the formalism to the 
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LES framework. To this aim, a detailed Lagrangian analysis of the error should 
be performed to complement present results and assist in model development. 
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Table [H Particle parameters; is the particle relaxation time, dp the particle 
diameter, Vg the particle settling velocity and Rcp the particle Reynolds num- 
ber. The superscript + denotes the corresponding quantities expressed in wall 
units. 
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St Tp (s) d+ dp ifim) V+ = g+ -St Re+ = V+ ■ d+/u+ 
1 1.133 • 10^3 0.153 20.4 0.0943 0.01443 
5 5.660 • 10^3 0.342 45.6 0.4717 0.16132 

25 28.32 • 10-3 0.765 102.0 2.3584 1.80418 

TABLE I. Particle parameters; Tp is the particle relaxation time, dp the particle diameter, 
Vs the particle settling velocity and Rep the particle Reynolds number. The superscript 
+ denotes the corresponding quantities expressed in wall units. 
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St = 5 particles with cut-off filter. Open symbols are used for the computed 
PDFs (o: CF=2, □: CF=4, V: CF=8); lines for the corresponding Gaussian 
PDFs ( : CF=2, : CF=4, ■ ■ ■■: CF=8). 

Figure (TO) Probability density functions of the wall-normal component of 
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(o: St = 1, □: St = 5, V: St = 25); lines for the corresponding Gaussian 

PDFs ( : St=l, : St = 5, ■ ■ ■■: St = 25). 

Figure [12) Probability density functions of the wall-normal component of the 
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FIG. 1. Mean values of the filtering error 6u in the streamwise (a), spanwise (b) and 
wall-normal (c) directions as a function of at varying particle inertia. Profiles 
refer to the cut-off filter with CF=4. 
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Streamwise direction, x 



FIG. 2. Fluid streaks in the near-wall region. Low-speed (resp. high-speed) streaks 
are rendered using colored contours of negative (resp. positive) streamwise fluctuat- 
ing velocity u'^ on a horizontal plane at z'^ = 10 from the wall. Panel (a) shows the 
streaky structures obtained from DNS; the instantaneous distribution of the St = 25 
particles comprised between z'^ = 10 and the wall is also shown to highlight prefer- 
ential distribution into the low-speed streaks. Panels (b) and (c) show the streaky 

structures obtained from a-priori DNS with coarsening-factor 4 and 8, respectively. 
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FIG. 3. Effects of filtering on the mean wall-normal fluid velocity seen by the 
particles at varying particle inertia. Profiles refer to the cut-off filter with CF=4. 
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FIG. 4. Rms of the filtering error 6u in the streamwise (a), spanwise (b) and wall- 
normal (c) directions as a function of at varying particle inertia. Profiles refer to 
the cut-off filter with CF=4. 
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FIG. 5. Skewness and flatness of tlie filtering error 5u in the streamwise (a-b), 
spanwise (c-d) and wall-normal (e-f) directions as a function of z'^ at varying particle 
inertia. Profiles refer to the cut-off filter with CF=4. 
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FIG. 6. Mean values of the filtering error 6u in the streamwise (a), and wall-normal 
(b) directions as a function of z'^ at varying filter type and width. Profiles refer to 
the St = 5 particles. 
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FIG. 7. Rms of the filtering error bw. in the streamwise (a), spanwise (b) and normal 
(c) directions as a function of z'^ at varying filter type and width. Profiles refer to 
the St = 5 particles. 
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FIG. 8. Skewness and flatness of the filtering error 5u in the streamwise (a-b), 
spanwise (c-d) and wall-normal (e-f) directions as a function of at varying filter 
type and width. Profiles refer to the St = 5 particles. 
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FIG. 9. Probability density functions of the streamwise component of the filtering 
error for different filter widths. Profiles refer to results obtained for the St = 5 
particles with cut-off filter. Open symbols are used for the computed PDFs (o: 

CF=2, □: CF=4, V: CF=8); lines for the corresponding Gaussian PDFs ( : 

CF=2, : CF=4, • • ••: CF=8). 
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FIG. 10. Probability density functions of the wall-normal component of the filtering 
error for different filter widths. Profiles refer to results obtained for the St = 5 
particles with cut-off filter. Open symbols are used for the computed PDFs (o: 

CF=2, □: CF=4, V: CF=8); lines for the corresponding Gaussian PDFs ( : 

CF=2, : CF=4, • • ••: CF=8). 
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FIG. 11. Probability density functions of the streamwise component of the filtering 
error for different particle inertia. Profiles refer to results obtained using cut-off 
filter with CF=4. Open symbols are used for the computed PDFs (o: St = 1, □: 

St = 5, V: St = 25); lines for the corresponding Gaussian PDFs ( : St = 1, 

: St = 5, • • ••: St = 25). 
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FIG. 12. Probability density functions of the wall-normal component of the filtering 
error for different particle inertia. Profiles refer to results obtained using cut-off 
filter with CF=4. Open symbols are used for the computed PDFs (o: St = 1, □: 

St = 5, V: St = 25); lines for the corresponding Gaussian PDFs ( : St = 1, 

: St = 5, • • ••: St = 25). 
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